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£S) ' Critical scattering analyses for dilute antiferromagnets are made difficult by the lack of predicted 

^ | theoretical line shapes beyond mean-field models. Nevertheless, with the use of some general scal- 

■ ing assumptions we have developed a procedure by which we can analyze the equilibrium critical 

scattering in these systems for H = 0, the random-exchange Ising model, and, more importantly, 
for H > 0, the random-field Ising model. Our new fitting approach, as opposed to the more con- 
ventional techniques, allows us to obtain the universal critical behavior exponents and amplitude 
ratios as well as the critical line shapes. We discuss the technique as applied to Feo.93Zno.07F2- 
The general technique, however, should be applicable to other problems where the scattering line 
<—j ■ shapes are not well understood but scaling is expected to hold. 

G : 

: 

d . I. INTRODUCTION 

£: 

Characterizations of the critical behavior of model systems through experiments and simulations are essential to 
verify the validity of theoretical models of phase transitions. Scattering techniques are invaluable for characterizing 
the staggered magnetization (the order parameter), M s , the antiferromagnetic fluctuation correlation length, £, and 
O [ staggered susceptibility, x S i a s a function of temperature in pure and dilute antiferromagnets, which prove to be ideal 
physical realizations of many model systems. Neutron scattering has been particularly instrumental in studies of 
I ' antiferromagnetsa, though magnetic x-ray scattering has also been employed to a limited extentH. Likewise, pulsed 
. heat and optical techniquesa have been essential in determining the critical behavior of the specific heat, C m . The 
C*") ' Ising model is the simplest of systems, with each spin having only two possible states, and becomes an exact model for 
the anisotropic antiferromagnets as the temperature T approaches the transition temperature T c . Three of the most 
fundamental phase transitions are the pure Ising model, the random-exchange Ising model (REIM) and the random- 
field Ising model (RFIM). The pure Ising model has the exact OnsagerQ solution for dimension d — 2. For d = 3 
only approximate renormalization solutions, simulations and expansions exist. Antiferromagnets exist with strong 
anisotropy, WC-h as FeFi for d = 3. This system exhibits universal Ising static critical behavior close to the transition 
temperatureHTj. The random-exchange model system is realized when the magnetic ions are randomly substituted 
by diamagnetic ions in Fe x Zri\- x Fi for x > x p , where x p — 0.246 is the percolation threshold for this body-centered 
tetragonal magnetic lattice if only the dominant next-nearest neighbor interaction!! is considered. Below x p no phase 
transition can occur for geometric connectivity reasons. The anisotropw-,increases as x decreasesO. As a result of 
the high anisotropy for these systems measured using neutron scatteringcSl, asymptotic Ising critical behaviot,is well 
followed in the reduced temperature range \t\ = |(T — T c )/T c \ < 10~ 2 , where T c is the transition temperaturea. 

Various experimental techniques cap_be used to extract universal Ising parameters associated with the asymptotic 
critical behaviors in antiferromagnetstil The universal parameters accessible through scattering techniques include 
the exponents and amplitude ratios associated with the asymptotic power law behaviors 
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z=zt\tr , (i) 

where + and — refer to t > and t < 0, respectively, 

M, = M \t\ f> , (2) 

where M Q is nonzero for t < only, 

x s = x±\tr , (3) 

and, for antiferromagnets with quenched randomness, the disconnected susceptibility 
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xf s =x? s± \tr (4) 

for q 0. Pulsed heat and optical techniques can be used to obtain the specific heat critical behavior, 

C = A ± \t\- a +B , (5) 

which becomes symmetric and logarithmic, 

C = A\n\t\ , (6) 

when a — > 0. 

Equations 1-6 represent the asymptotic behaviors exhibited for data with sufficiently small \t\. The range of 
asymptotic Ising behavior in anisotropic short-range interaction antiferromagnets, for example, is partly determined 
by the anisotrop^-sisength. For example, whereas FeFi showsB asymptotic behavior for \t\ < 10~ 2 , the less anisotropic 
isomorph Mni^liirEj shows asymptotic behavior only for \t\ < 1CP 3 . If data are taken outside the asymptotic range of 
\t\, but still within the critical region, fits to the power law expressions yield only effective exponentsEfl and effective 
amplitude ratios. In such cases, it may be more effective to use scaling function analyses, which include crossover to 
asymptotic critical behaviors, when some quantity such as the applied field, H, can be varied. We will briefly discuss 
this with regard to specific heat analyses done principally using Feo. 93.Zno.07F2 data. Scaling functions are also of 
great utility in fitting scattering data when the line shapes are not well known. For example, at constant H in the 
asymptotic region, we can use the fact that the scaling functions can depend only on \q\/ k., where q is the distance 
from the antiferromagnetic Bragg scattering point, (100) for FeF2, in reciprocal lattice units (rlu) and k = l/£ is the 
reciprocal correlation length for antiferromagnetic fluctuations. In this work, we will show how the general properties 
of the scaling functions caa be utilized to advantage in the characterization of the random-exchange and random-field 
scattering data obtainedcJ using Feo. 93-^^0.07^2- However, the technique has more general utility; it can be used in 
any case where theory does not provide adequate models of scattering functions but scaling is expected to hold. 



II. SPECIFIC HEAT SCALING 



Although we primarily focus here on the scaling behavior of the critical scattering of the d = 3 dilute antiferromag- 
net, it is instructive to review briefly the success of scaling analyses in the study of the specific heat critical behavior in 
the same system. Not only will this demonstrate the usefulness of the scaling approach, it will highlight the results for 
the specific heat behavior, which are complementary to the scattering results but show the largest discrepancy with 
simulation results for the RFIM. This will be important in our later assessment of the agreement between simulations 
and experimental results. 

In zero field, the dilute anisotropic antiferromagnet is predicted to have a transition described by the random- 
exchange Ising model. One of the most striking changes in the critical behavior induced by random quenched dilution 
is observed in the specific heat for d = 3. Whereas the pure FeF% sample shows a = 0.11, in agreement with theoryli3, 
the Harris criterion imposes the constraint that a < upon dilution. Indeed, tiie exponent has been foundlill to be 
a = —0.10 ±0.02 experimentally. Monte Carlo studies also yield a negative valueH. Interestingly, with the application 
of a field, the RFIM specific heat is again found to be divergent in experimentsEZl, with a » 0, as discussed below. 

One way to utilize scaling functions is to attempt to collapse the experimental data onto a scaling function of the 
appropriate scaling argument. The collapse will only work well if the critical parameters used in tiie data collapse are 
correct. The RFIM scaling behavior of the free energy, for example, is expected to have the formO 

F~ff 2 ( 2 - Q )/^ 5 (|t H |ff- 2 /^) , (7) 

where tjj = (T — Tjv + bH 2 )/TN, Tn is the zero-field transition temperature, b is a small mean-field (MF) parameter, 
a is the zero- field specific heat exponent, and 4>rf is the RFIM crossover exponent. When a new phase transition 
occurs at T = T C (H), the asymptotic limit of the specific heat can be obtainecO from the free energy in the limit 
1*1-0 

B 2 F 

C = — = H- 2a ^ RF g'(\t H \H- 2 ^ RF ) ~ #»(a-°0/W|t|-a , (8) 

where t = (T - T C {H))/T C (H). For a = 0, this becomes 

C - H- 2a/<t,RF In \t\ , (9) 
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which is symmetric above and below T C (H). Note that the field dependence of the peak amplitude is dependent on 
a. In the scaling plots, both the shape and the field dependent amplitude must be correct for the data to collapse 
onto a single scaling function. 

The scaling behavior of the specific heat has been experimentally demonstrated for both d — 2, where no transition 
takes placeE3, and d — 3, where a new transition occursliH. In both cases, the critical parameters can be determined 
accurately by the quality of the data collapse onto a single scaling function. The most accurate measurements 
have been obtained using the octical linear birefringence technique which has been shown to faithfully represent the 
magnetic specific heat behaviomO'EiJ. It was shown that the experimental data for the d = 3 system Feo.g3Zno.07F2, 
when divided by u~ 2a ^ RF , collapse onto a single scaling function if a = — 0.10±0.02 and <Prf = J-|42±0.03, the laiter 
having been determined by measurements on Fe x Zni_ x F 2 with several different concentrationsE3 and predicted^ to 
be a few percent larger than the zero-field staggered susceptibility exponent measuredE-3 to be 7 = 1.31 ± 0.03. Note 
that if the RFIM specific heat exponent a — 0, as indicated by the asymptotic shape of the curve, the random- 
exchange-exponent must have the value a = —0.10 ± 0.02, which is consistent with the earlier value a = —0.09 ± 0.03 
obtainedc3 by fitting the data from a lower magnetic concentration sample with x = 0.60 to a power law. Hence, both 
the shape and field-dependent amplitude of the critical peak are consistent with a = 0, i.e. a symmetric, logarithmic 
divergence. 

Further evidence that a is close to zero is obtained from Faraday rotation experiments under constant field or 
constant temperature. Faraday measurements yield the critical behavior of the specific heat, but with_a different field 
amplitude than the specific heat peaks. Near the phase transition it was shown by Kleemann, et alH3, that 
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which becomes, for a = 0, 



jj2(l+a~a-4> RF /2)/4> R F^-a /-^q\ 



and 



_ d 2 F Tr2(2+a-a-4>iu*)/<t>RFU\-a fi o\ 

dH ~ 3H 2 m ' [ ' 

which becomes, for a = 0, 

^„fT2(2-"-*KF)/***]n|t| . (13) 

The field dependent amplitudes of logarithmic peaks have been measured with the results 2(1 — a — 4>rf /2)/4>rf ~ 
0.56 and 2(2 — a — (f>RF)/4>RF ~ 0.97. With the measured value <Prf — 1-42 ± 0.03, these two equations yield 
a = —0.11 ± 0.02 and a = —0.11 ± CU34, respectively. These values are consistent with the values_p = —0.10 ± 0.02 
from the Feo.g3Zno.07F2 experimentll3 and a = —0.09 ±0.03 from the Feo.eoZno.40F2 experimented. Hence, the field 
dependent amplitudes as well as the peak shapes are all consistent with a « 0. A scaling analysis would similarly 
require a « in the case,of Faraday rotation for a good data collapse. The result a « from experiments contradicts 
the Monte Carlo result^ a = -0.5 ± 0.2. 



III. SCATTERING SCALING FUNCTION 



We next turn to the scattering function, which should obey scaling properties close to T C (H). Within the static 
approximation^!, the intensity of the magnetic critical scattering from high quality single crystal magnetic systems is 
proportional to the Fourier transform of the spin-spin correlation function S(q) = [< s q s- q >] convoluted with the 
instrumental resolution, where the <> brackets signify a thermal average and the [ ] brackets signify a configurational 
average. 

The spin-spin correlation can be expressed as 

S(q) = X s + xf S , (14) 

where 

X s = [< SqS- q > - < S q >< S_, >] (15) 
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is the staggered susceptibility and 

X f s = [< s q X s. q >] (16) 

is the disconnected susceptibility. For line shapes obtained at one value of the field, we expect a scaling function that 
only depends on the ratio \q\/ n of the two physically relevant inverse length scales. For |g| > 0, 

X.{q) = A ± K"- 2 f{q/K) . (17) 

For pure, translationally invariant systems xf ls — Ms 2 S(q). For random systems, on the other hand, xt ls mayJiave 
a q dependent contribution. In particular, in RFIM systems, such a term is induced by the random field, givingtil for 
| q | > at constant H . 

xf s (q)=A ± B ± K f )- i g(q/ K ) . (18) 

Hence, S(q) involves two possibly independent scaling functions, f(q/K) and g(q/K) and two possibly independent 
exponents 7 and 7. This makes the method of collapsing data onto scaling functions extremely difficult since the 
contributions to the data from the two scaling functions are not easily separable. 

To fit the data we must, in principle, use data only in the range of small k (i.e. small \t\) and \q\ in order to 
be sure we are dealing with asymptotic behavior. In the study of FeF2, data were used in the range \t\ < 10 -2 to 
obtain the most reliable exponents and amplitude ratios, since this was shown to be in the asymptotic range for pure 
Ising behavior in specific heat critical behaviorU measurements. In the Feo.g3Zno.orF2 experiments, we used data 
for \t\ < 1CP 2 since this is the range for which the specific heat shows the RFIM logarithmic behaviorEZl. However, 
we typically use a wide range of \q\. The crossover at large \q\ is relatively unimportant since the critical scattering 
intensity becomes very small. On the other hand, including data at large \q\ helps to set the level of background 
scattering. 

Many scattering critical behavior analyses are done using the simple MF Lorentzian for \q\ > 0, 

f(q) = , A± 2/ 2 > ( 19 ) 

1 + q I K 

giving 

S(q) = Xs (q) = -^—2 (20) 
q z + K 

for q ^ 0, consistent with the fact that 77 = in MF. Since the upper critical dimension, d u , above which the MF 
equations are correct, is four or greater, the Lorentzian line shape can only be approximate for three dimensions. 
Deviations from the MF Lorentzian line shape should become evident as one approaches the transition temperature 
and are generally found to be much more important .below the transition temperature. This has been discussed with 
respect to the pure d = 3 Ising antiferroniagnet FeF2^. Scattering data for FeF2 were analyzed using approximations 
to the line shapes by Fisher and BurfordElT(BF) for T > T N , 

(l + ^y/^y/ 2 

fiq/K) i + ^ 2 /. 2 ' (21) 

and by Tarko and FisherS (TF) for T <T N , 

where 0, 0', 0", a, tp = 1 + l/2?70 2 and V' = 1 + l/2#' 2 + er(0' 2 - ft' 2 ) are fixed to values determined from the 
numerical studies. The values are given in Table I. The expressions have the correct scaling behavior in the limits 
\q\/K — > and \q\/n — > 00 and serve as appropriate interpolation functions between those limits. We show the critical 
exponents and amplitude ratios obtained from fits of the data obtained for 10 -4 < \t\ < 1CP 2 using these scaling forms 
in Table I as well as theoretically determined universal critical parameters. The experimental and theoretical values 
serve to contrast those from similar analyses done on the diluted system Feo.93Zno.07F2 to be discussed next. The 
corrections to MF are only significant for \t\ < 1CP 3 and are most significant for T < TV. Since many d = 3 studies 
do not probe critical behavior any closer than this to T c , the Lorentzian line shape generally serves satisfactorily for 
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extracting estimates of the critical exponents. Obviously, more precise measurements that probe regions of smaller 
\t\ can yield much more accurate critical parameters, but only if a suitable line shape, is used. 

Data for the pure d = 2 case of K2C0F4 has been analyzed in a similar mannerEa. In this case rj = 0.25 and an 
analysis using the MF Lorentzian fails markedly below Tjy. The appropriate TF and FB equations, in contrast, yield 
critical behavior consistent with the d = 2 Ising model. 

The scattering function for the dilute antiferromagnet is predicted to have an additional term not present in the 
pure caseE3 that may affect experimentally determined to corrections to scaling and, possibly, measurements of the 
amplitude ratio for \ s - The exponents are probably not influenced since this extra contribution vanishes as T — > T/v. 
A fit to a simple Lorentzian seems to work reasonably well for data with |f | < 10~ 3 in the very dilute antiferromagnet 
Feo.4sZno. 54^2- For the random-exchange Ising modeO, r\ « 0.04, which is similar to the value for the pure case. 
Hence, it is reasonable that the Lorentzian line shape works well in this reduced temperature range just as it does in 
the pure d = 3 Ising case. For data closer to T c , however, just as in the pure case, we would expect deviations from 
the MF Lorentzian, particularly for T < T c . Unlike the Tarko-Fisher case for the pure Ising model, no approximants 
have been worked out for the random-exchange model for use in data analyses beyond MF. One possible approach 
is to use the same TF and FB expressions developed for the pure case. This assumes that the line shapes for the 
pure and REIM are very similar. Such an approach was successfully employed to analyze the data for the dilute 
d — 2 antiferromagnetEi Another strategy is to use the same forms, but to let the TF/FB parameters be free fitting 
parameters. Such functions would satisfy scaling in the proper limits and would hopefully be very good interpolative 
functions between those limits. We will describe below results for Feo. 93^^0.07-^2 obtained using the latter of these 
two methods. 

We next turn to the more difficult case of the random-field scattering in Feo^Zno Q7F2, which occurs for H > 0. 
In the case studied H = IT . We expect the scattering function for |g| > to be 

S(q) = A ± K r '- 2 f{q/K) + B ± A ±2 K ^- 4 g{q/ K ) . (23) 

Taking into consideration the instrumental resolution as well as the two separate scaling functions, there is little 
chance of using the data directly to determine the two independent scaling functions. Hence, we must start with 
model functions_and test their appropriateness. The first natural test functions to use in the data analysis are the 
MF scaling oneso, 

S ( q > = 1l2 /2 + ( 1 i 2/ 2\2 ' ( 24 ) 
1 + q z I n z {\ + q z lK z Y 

where we might expect the amplitudes A^ and B^ to be temperature dependent, since rj and fj are not, in fact, zero. 
We did such an analysis previouslylEl and found reasonable fits at all temperatures in the sense that the line shapes 
yielded values of k and x s = A^k^ 2 . However, when an attempt was made to fit these values to power law behaviors, 
reasonable results were obtained above the transition but not below. The results above T C (H) are v — 0.90 ±0.01 and 
7 = 1.72 ± 0.02. These values are consistent with results above T C (H), obtaineda for x — 0.6 using a MF analysis. 
At these lower concentrations, equilibrium critical behavior is not obtained below T C (H), but the data analysis was 
done well above T C (H) where equilibrium prevails. I— ■ 

We note that the failure of the MF analysis below T C (H) is very similar to the situation observed^ in the pure 
d = 2 antiferromagnet K2C0F4. The MF equations fail in that case because d = 2 is far from the upper critical 
dimension d u — 4 and r\ — 0.25 in contrast with the pure d = 3 case where 77 = 0.04. The TF expression developed for 
d = 2 served nicely for the data analysis and agreement was found with theory. Interestingly, the FB scaling function 
for T > Tm is not very different from the MF one. For the d = 3 RFIM, the value of r\ is predicted to be even larger, 
perhaps as large as r\ = 0.5. Hence, it is not that surprising that the MF data analysis fails below T C (H) in this case. 
To proceed, we must go beyond a simple MF line shape analysis. 

Since there is not yet a theoretical line shape available beyond MF, we must try to use scaling properties to guide 
us. However, if there are two independent scaling functions, as in Eq. 24 r .the task becomes formidable. Fortunately, 
there are two approximations motivated by theoreticaO and simulatioriifl works. The first is that fj — 2rj, a limiting 
case of the Schwartz-SoffeiEj inequality fj < 2t]. The second is that g(q/n) = f 2 (q/ k) is a very good approximation. 
We adopt these two simplifications, making the scattering function for \q\ > the more manageable expression 

S(q)=A ± K'l- 2 f(q/K)(l + B ± A ± ^- 2 f(q/ K )) . (25) 

It is not clear at this point how accurate these approximations are. However, it is not possible to proceed without 
them and they appear to be well justified. It is highly unlikely that experiments will be able to test the validity of 
the assumptions directly. Only further theoretical progress can provide a better starting point for the data analysis. 
Note that the MF expression in Eq. 24 is a special case of this for 77 = 0. We are still faced with the correction to 
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instrumental resolution, however, which itself depends on the line shape we are trying to determines. Hence, it is 
still hard to scale the data directly without an explicit functional form for f(q/n). To simplify the procedure further, 
therefore, we have adopted as our scaling functions the TF and FB scaling functions except that we allow all the 
parameters to vary. Hence, we are assured of the proper scaling in the limiting cases of \q\ — > and k — > and 
hopefully the interpolation between these limits will be adequate with the parameters determined from the fits of the 
data. Unfortunately, the usual technique of fitting each scan separately, where a scan is made in q at fixed T, and 
then extracting the exponents from power law fits to the resulting k and \ s cannot work well since there are now so 
many free parameters, including the various exponents, amplitudes and the TF/FB parameters. We can, however, fit 
all the data scans simultaneously, since the line shape parameters are all the same for every temperature for T < T c 
and T > T c and the critical exponents are the same for all T. The new technique has the advantage over the more 
classic technique in that the line shape does not need to be known beforehand. In the RFIM, the line shape is both 
unknown and far from the MF prediction, so the classic method using the MF line shapes failed to yield the critical 
behavior parameters. Our new procedure using the TF and FB line shapes with variable parameters and fitting all 
the data simultaneously, on the other hand, was successfully employed for the H > random-field Ising behaviors in 
Feo.93Zno.07F2 and yielded both the critical behavior parameters and the critical scattering line shapes. Both the 
classic and new techniques worked well for the H = case in which the line shapes are nearly mean-field and much 
simpler than the RFIM ones. This demonstrates the reliability of the new technique and we can apply it with some 
confidence to the RFIM, where the more classical technique fails to yield results. 



IV. EXPERIMENTAL AND FITTING DETAILS 



Adding to the difficulty of implementing scaling in the scattering line shape analysis is the necessity to account 
for significant instrumental resolution corrections. The instrumental resolution can be measured for a particular 
spectrometer configuration by measuring the width of the magnetic Bragg scattering peak at low temperatures. The 
Bragg scans in the transverse, longitudinal and vertical directions well below the transition temperature yield the 
response to scattering from the Bragg peak. Theoretically, the Bragg peak is a delta function and in practice is 
much more narrow than the instrumental resolution in good quality crystals. Since the scattering in these anisotropic 
crystals is well described within the static approximation, the scans in the three directions yield the widths along the 
three principal axes of the resolution ellipsoid. If we have a theoretical line shape, we can use the measured Bragg 
scattering scans to numerically integrate the line shape which can then be compared directly with the scattering data. 
Alternatively, other groups have used Gaussian (or with less accuracy triangular) approximations to the Bragg scans 
and then analytically integrated to obtain the resolution corrections. In this study we exclusively use the numerical 
integration technique. For resolution curves measured at uniform steps in q, we have for the intensity in transverse 
data scans, 

J( , J2 S((g- go- a) 2 + b 2 +c 2 )T a L b V c 

m ET^v c ' (26) 

where the sums are over a, b, and c and qo accounts for imperfect alignment in the transverse direction. For a well 
aligned crystal, qo is usually much smaller than the resolution width in the transverse direction. Misalignments along 
the vertical and longitudinal directions are generally inconsequential for a well aligned crystal since the resolutions 
are much courser in these directions than in the transverse one. T a , Lb and V c are the approximately Gaussian Bragg 
line shapes measured at low T using evenly spaced steps in q. This is the technique explained!! in detail in the context 
of the study of FeF^.. If a crystal has a small mosaic, this can also be approximately taken into account in the same 
manner. 

To fit the resolution-convoluted line shape to the scattering data, we use a nonlinear least squares fitting routine 
to determine the parameters of the theoretical line shape and critical exponents. Note that with each iteration of 
the fitting program, the line shape must be reintegrated over the resolution ellipsoid since the amount of correction 
from the resolution convolution depends on the line shape and, hence, the fitted parameters. Since the resolution 
correction in neutron scattering experiments is substantial and line shape dependent, it is very difficult to determine 
the line shape directly from the data. Hence, in the absence of a theoretical model, we must choose a trial function 
that satisfies the correct scaling requirements and allows suitable flexibility in the fits of the data. Our strategy is to 
use the TF and FB line shapes described above. 

For the ex-oeriments we used two samples of Feo.93Zno.07F2- One is the same large sample used in specific heat 
experimental^. It is somewhat irregular in shape and has a mass of 1.35 g. The magnetic concentration gradient 
limited the range of data unaffected by rounding to \t\ > 1.15 x 1CP 3 . The second is a slice cut from the large sample 
with its faces perpendicular to the magnetic concentration gradient. It is approximately one tenth the mass of the 
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original sample. The smaller sample was used to obtain data closer to the transition with \t\ > 1.14 x 1CP 4 . The 
data from both samples were used simultaneously in the data fits with different instrumental resolution corrections 
appropriate to the two samples and the spectrometer configurations used to make measurements on them. Neutron 
scattering measurements were made at the Oak Ridge National Laboratory High Flux Isotope Reactor using a two- 
axis spectrometer configuration. We used the (0 2) reflection of pyrolytic graphite (PG) at an energy of 14.7 meV 
to monochromate the beam. We mainly employed two different collimation configurations. The lower resolution, 
primarily used for the large sample, is with 70 min of arc before the monochromator, 20 before the sample and 20 
after the sample. Primarily for the thin sample, we made scans with 10 min of arc before and after the sample. Two 
PG filters were used to eliminate higher-order scattering. The carbon thermometry scale was calibrated to agree with 
recent specific heat resultstZl for the H = transition. The field dependence of the thermometry was also calibrated. 
All scans used in this report, other than those used for obtaining the resolution ellipsoid, are transverse ones about 
the (10 0) antiferromagnetic Bragg point. 

V. FITTING RESULTS 

We applied the techniques described above to the random-exchange behavior (H = 0), where there is only one term 
in the scattering function for |g| > 0. We first did the fits using the MF line shapes imposed by setting a = 1 and 
(j), <j)' and <f>" equal to zero. The results for the fitted critical parameters arc essentially identical to the results from 
the conventional method of fitting each scan to obtain the temperature dependent correlation and susceptibility and 
subsequently fitting them to extract the exponents and amplitude ratios. The results are shown in Table I. We include 
the expressions for the normalized x 2 m the table, which are most useful for relative comparisons between fits. The 
values are not close enough to unity to consider the fitting functions statistically perfect ones. Very small systematic 
errors from the resolution corrections and the approximations in the line shapes can easily account for values of x 2 
being somewhat larger than unity. The goodness of the fits are better judged from the scaling plots discussed below. 
The agreement between the two methods of fitting the data gives us some confidence in fitting the field data, for 
which the line shape is unknown, using the technique in which all data are simultaneously fit. We proceed beyond 
MF by doing the fits with the TF and BF parameters as free parameters. The fitting results are again shown in 
Table II. The range of data was restricted to the reduced temperature range \t\ < 0.15 and includes 2198 data. The 
results for the critical exponents and amplitudes are not very different from the previous result and the line shape is 
not very different from the MF one. When examining individual scans, we observe no systematic deviations of the 
data from the fit at any temperature for which data were included in the fit. We show representative scans in Fig. 
1, including data above and below T c . No systematic deviations of the data from the fits are evident, indicating that 
the line shapes work well for the data within the range of \t\ used in the fits. To demonstrate that the data are well 
described by the fitted scaling function, we subtract the scattering background from them, deconvolute them with 
the instrumental resolution, divide them by A^n 11-2 and plot them versus \q\/k in Fig. 2. The data from the two 
samples are plotted separately, but the solid curves are identical in the two cases. The scatter of the data for the 
small sample simply reflects the smaller size and the resulting lower count rate. However, some of the data from the 
small sample are taken much closer to T c , since the rounding from the concentration gradient is less significant. The 
data close to T c are very important in the fits. The consistent results obtained using MF and scaling line shapes in 
the analyses give us confidence that the scaling technique may work well even in cases where the line shapes are very 
far from being MF, as in the RFIM case we describe next. 

For the RFIM case, the procedure is identical to the random-exchange one above, except that the scattering function 
involves the more complicated expression in Eq. 25. In one fit, we restricted the data range to \t\ < 10~ 2 , with 2444 
data, since this is the temperature range over which the asymptotic logarithmic specific heat behavior is observedO. 
In a second fit, we further restricted the temperature range to \t\ < 3 x 10~ 3 , with 1000 data, to test whether 
crossover effects were still significant. The fitting results are shown in Table III. The two fits are rather consistent, 
suggesting that the results are close to the asymptotic ones. In Fig. 3 we show data and the corresponding fits for a 
few scans to show that there are no significant systematic deviations. After determining the line shape parameters 
and critical exponents, we used a procedure similar to that described above to demonstrate that the data collapse 
onto scaling functions for T < T c and T > T c . We subtracted the scattering background, deconvoluted the data with 
the instrumental resolutions, divided by v4 ± k 1, ~ 2 (1 4- B ± A ± K v ~ 2 f(q/K)) and, for clarity, plotted the results versus 
K ± |g|/K instead of simply \q\/n so that the data for T > T C (H) and T < T C {H) do not overlap. The results of this 
procedure are shown in Fig. 4. The data from the large and small samples are shown separately. 

We have obtained from the scaling analysis not only the critical parameters but adequate approximations to the 
line shape scaling functions for both the d — 3 random-exchange and random-field Ising models. We already have 
the d = 2 and d — 3 pure line shapes from the TF and FB expressions. It is instructive to compare all of these to 
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the MF Lorentzian line shapes. The comparisons are shown in Fig. 5, where we have plotted scaling functions versus 
The upper plot is for T > T c and the lower one for T > T c . The lowest curve in each case represents the 
simple Lorentzian line shape that is accurate in the MF limit. The pure d — 3 Ising model line shapes show nearly 
MF behavior for T > T c . For T < T c , more significant deviations from MF are apparent. This is consistent with 
the experimental results obtained using FeF 2 . The d = 3 experimental random-exchange scaling functions indicate 
that there is not much difference between the pure and random-exchange line shapes for d = 3. In contrast, for the 
random-field case, the deviations from MF found in the experiments are very large for T < T C (H) in comparison to 
the pure and random-exchange d = 3 Ising models. Even in the case for T > T C (H) the deviations are large relative to 
the other d = 3 models. The very large deviations for T < T C (H) are consistent with the failure of the MF Lorentzian 
line shape in producing values for k and \ s that obey power law behavior. For the pure d — 2 Ising model, shown as a 
reference for non-MF behavior, the results for T < T C (H) are the largest of the examples shown, whereas the behavior 
for T > T c is not far from MF. We see from this comparison of the scaling functions, particularly for T < T c , that MF 
behavior is fairly well followed for the d = 3 pure and random-exchange models, whereas for the random-field case, 
where the upper critical dimension has significantly increased to d u — 6, the line shapes are very non-Lorentzian. 
Determining the line shape from the scattering data is difficult. We have achieved an approximate determination of 
the proper line shapes. Once theoretical results give a more firm foundation for the trial function for the scattering 
line shape, we will be able to give a more concrete comparison of the experiments to theory and simulations. 

In Table IV we show the results of the critical exponents for the specific heat and neutron scattering RFIM 
experiments compared with simulation results. For the pure d = 3 Ising model, the experimentally determined 
universal exponents and amplitude ratios agree very well with theoretical and simulation results. The experimental 
results from neutron scattering were obtained using the TF and FB line shapes. The specific heat was obtained using 
both pulsed heat and birefringence techniques and clearly can be considered to be exceedingly well characterized. For 
the REIM, we again have excellent agreement between experiment and Monte Carlo simulations. The specific heat is 
determined most precisely using the birefringence technique, but the results are consistent with the pulsed heat data. 
This case can also be considered well characterized. 

The RFIM case shows mixed results when the experimental exponents are compared with those from simulations. 
The scattering results for v and 7 are quite consistent with the simulation values. The exponent has not been 
reliably measured yet and its comparison with simulation results is very important. The most glaring inconsistency 
is between the simulation and experimental values of a, where the experiments indicate a symmetric, logarithmic 
divergence and the simulations indicate a non-divergent peak. This inconsistency deserves further study. 

Although we have achieved fits to the data obtained in equilibrium scattering experiments using Feo.g3Zno.01F2, 
the asymptotic behavior is only observed very close to T C (H) at H = 7 T. Our results may be somewhat influenced by 
the effects of crossover to random-exchange behavior. We hope in the future to do these experiments at much higher 
fields to ensure that the results we have obtained in this study are close to the asymptotic ones. In addition, by using 
apparatus capable of reaching much higher fields, we can study the scaling behavior of the scattering line shapes as 
a function of \t\H~ 2 ^ RF in a similar way to our treatment of the specific heat. 

This work has been supported by DOE Grant No. DE-FG03-87ER45324 and by the Oak Ridge National Laboratory, 
which is managed by UT-Battelle, LLC, for the U.S. Dept. of Energy under contract DE-AC05-00OR22725. 
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-0.1 -0.05 0.05 0.1 
q (rlu) 

FIG. 1. The logarithm of the scattering intensity vs. q and the curves representing the corresponding fits for H = 0. One 
scan for T > T C (H) and one for T < T C (H) are shown for the large sample in figure a and similarly for the small sample in 
figure b. 
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FIG. 2. Scaled neutron scattering data, deconvoluted with the instrumental resolution, taken at different temperatures at 
H = T collapsed onto the universal function /(<?/«). The scatter in the small sample data is larger due to smaller number 
of counts obtained in the thin sample. The fit was made for |£| < 0.01. The solid curves, which are identical in the upper and 
lower plots, represent the line shapes for T > T c and T < T c determined from the fits to the data. 
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-0.1 -0.05 0.05 0.1 
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FIG. 3. The same as in Fig. 1, but for the the random-field Ising case H = 7 T. 
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FIG. 4. The same as in Fig. 2, but for the random-field Ising case H = 7 T. The fit was made for \t\ < 0.01. 
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FIG. 5. A comparison of the logarithm of the scaling functions versus q/n for different models. The pure cases are from the 
TF and FB expressions. The REIM and RFIM are from the experiments, as described in the text. Note that the corrections 
to the MF equation are largest below the transition. The random-field deviations for d = 3 are greater than the pure and 
random-exchange, but are significantly smaller than the pure d = 2 scaling function below the transition. Above the transition, 
all of the line shapes are close to MF. 



14 



w 



parameter FB and TF d = 3 Ising 

T c 78.418 ± 0.001 K 

v 0.64 ± 0.01 0.63002 ± 00023 

k+ /kq 0.53 ± 0.01 0.510 ± 0.002 

7 1.25 ±0.02 1.2371 ±0.004 

Xt/Xo 4.6 ±0.2 4.77 ±0.02 

t] 0.056* 

i> i.ooi* 

cr 2?7 = 0.111* 

0.15* 

</>' 0.3247* 

4>" 0.09355* 

4>' 1.0137* 



TABLE I. Experimental and Theoretical values for the critical parameters of the d = 3 pure Ising model. Data were fit 
over the range 10" 4 < \t\ < 10" 2 . Note that ip = 1 + l/2r)<f) 2 and tp' = 1 ± l/2ry0' 2 ± a((p' 2 - 0"' 2 ). Values marked with * are 
fixed parameters in the fitting process. 



parameter 


Lorentzian 


TF/FB 


T c (fixed) 


72.73iT 


72.737T 


V 


0.079 ±0.010 


0.079 ±0.012 


V 


0.70 ±0.03 


0.70 ±0.02 


A+ 


7.66 ±0.15 


7.59 ±0.10 


A~ 


6.71 ±0.10 


6.42 ±0.10 




0.57 ±0.02 


0.56 ±0.02 


K 


1.21 ±0.10 


1.13 ±0.05 


cr 


1* 


0.16 ±0.20 


<t> 


0" 


0.18 ±0.02 


4>' 


0* 


0.18 ±0.02 


4>" 


0* 


0.08 ±0.10 


c 


0.0142 ± 0.0002 


0.0142 ± 0.0001 




2.0 


1.7 


No. pts. 


2280 


2198 



TABLE 11. The values found for the parameters from the fits for the H — T data. The exponents and amplitude ratios 
are defined in the text. C refers to a constant, q-independent background scattering term. The scattering data for the large 
sample were fit for |t| > 1.15 x 10 -3 and the data for the small sample were fit for |t| > 1.14 x 10~ 4 . T c was fixed in the case 
of the scattering fits. The first column of results is for the MF line shapes. The second is TF/FB using the pure line shape 
parameters. The third column is obtained by letting the TF/FB parameters be fit along with the exponents. 
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parameter 


|t| < 10 


1 + 1 ^ O w 1 A - 3 

|t| < o X 1U 


i c (faxed) 


70.61a 


70.61A 


V 


0.20 ± 0.05 


0.16 ± 0.06 


V 


0.88 ± 0.05 


0.87 ± 0.07 


A 


i a a _i_ n o 


A 01 _l_ A O 


A 


6.15 ± 0.14 


A ATL 1 A 1 C 

4.45 ± 0.15 


+ 

Kg 


1.13 ± 0.04 


A A C I f\ 1*7 

0.95 ± 0.17 


Kg 


nod i /a -i -i 

3.24 ±0.11 


2.78 ± 0.5 




(AT _l_ A 1 \ w 1 A — 5 

(4.7 ± 0.1) X 10 


/ O AA 1 A1Q\w1A — 5 

(3.00 ± 0.13) x 10 


_fc> 


! A A 1 n o\ w 1 A - 5 

(4.0 ± 0.3) x 10 


/O A _l_ 1 A\ w 1 A — 5 

(8.0 ± 1.0) x 10 


O 


0.0 / ± 0.0 


0.80 ± 0.0 


<t> 


0.16 ±0.04 


0.08 ±0.01 


0' 


0.39 ±0.25 


0.36 ±0.3 




0.31 ±0.25 


0.26 ±0.2 


c 


0.017 ±0.001 


0.016 ±0.001 


X 2 


3.07 


2.3 


No. pts. 


2444 


1000 


TABLE III. 
were fit for \t\ 


The values found for the parameters from the TF/FB fits for the H — 7 T data. 
> 1.15 x 10~ 3 and the data for the small sample were fit for \t\ > 1.14 x 10~ 4 . 


The data for the large sample 
T c was fixed in the case of the 



scattering data. 



Pure 



7 



Random 

Exchange 

a 


v 

7 
'/ 



Random 
Field 



n 

7 



FeF 2 
0.11 ±0.00 M 
0.325 ± 0.0050 
0.64±0.0li 
1.25 ±0.02i 

o.osi 



Fe x Zm- x F2 

(H = 0) 
-0.10 ±0.0 
0.350 ± O.Oi 
0.69 ±0.0: 
1.31 ±(U13Efl 

o.ioy 



Fe x Zm- x F2 

(#>0Y-, 
0.0 ± 0.02O 

not measuredB 

0.88 ±0.05 

1.58 ±0.13 

0.20 ±0.05 
2 7 = 3.16 ±0.26 

2?7 = 0.40 ±0.10 



Renormalization and 
High T Expansie«|S 

0.1099 ±o.oooaia 

0.109 ±0.0040,-, 
0.32648 ± 0.000 

0.3258 ± 0.0014E3-, 
0.63002 ± 0.0002313 
0.6304 ±0.0013 21 
1.2371 ±0.0004 i£ 
1.2396 ±0.001322 
0.0364 ± 0.0004 i£ 
0.0335 ± 0.0025E! 



Monte Carlo 



-0.051 ±0.013 
0.3546 ± 0.0028^ 
0.6837 ± 0.005,3 

1.342 ±0.01C 
0.0374 ± 0.0045 



Monte Carlo 
Exact GrounidnState 
-0.5±0.2Ei 
-0.55±0.2S 
0.00 ±0.05^ 
0.02±0.OJip 

1.1 ± 0.2Ei 
1.14±0.ldp 

1.7±0.2q 

1.5 ± 0.2EjL 

0.50 ± o.oia 
3.3±o.eq 

3.4±0.4EjL 
1.03 ±0.05^1 



TABLE IV. The d = 3 Ising critical exponents obtained from experiments on the Fe x Zn\ — xFi and the corresponding 
values from Monte Carlo simulations and exact grounA state-ralculations. More extensive experimental, tjheoretiefil and simu- 
lation results are compared for the pure model in ref.EJ andE3 and for the random-exchange model iip-ref.H andcil. Note that, 
unlike the pure and REIM cases, it is well established that hyperscaling is violated in the RFIM casetil, i.e. a + dv 7^ 2 for the 
RFIM. 
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